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Abstract 



In this work we study the nonnegative solutions of the elliptic system 



Av = |x| V 



in the superlinear case ^5 > 1, which blow up near the boundary of a domain of R^, or 
at one isolated point. In the radial case we give the precise behavior of the large solutions 
near the boundary in any dimension N. We also show the existence of infinitely many 
solutions blowing up at 0. Furthermore, we show that there exists a global positive solution 
in R^\{0}, large at 0, and we describe its behavior. We apply the results to the sign 
changing solutions of the biharmonic equation 



Our results are based on a new dynamical approach of the radial system by means of a 
quadratic system of order 4, introduced in [4], combined with the nonradial upper estimates 
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u = \x\ \u\ . 



of O. 



1 Introduction 

This article is concerned with the nonnegative large solutions of the elliptic system 

An = \x\°-v^ I . 

A?; = |x|V, ^ ' 

in two cases: solutions in a bounded domain in M^, which blow up at the boundary, that is 

lim uix) = lim vix) = oo, (1-2) 

where d{x,dO,) is the distance from x to dQ; or solutions in Q\ {0} which blow up at : 

lim u{x) = oo or lim v{x) = oo. (1-3) 

We study the super linear case, where fi,5 > 0, and 

D = n6-1>0, (1.4) 

and a, b are real numbers such that 

a,5 > max{-2, -A^}. (1.5) 

First we recall some well-known results in the scalar case of the Emden-Fowler equation 

AU = U^ (1.6) 

with Q > 1. Concerning the boundary blow-up problem, there exists a unique solution U in Q 
such that lim^(2, gf2)__^o ^(^) = i^^ar 

?7(x) = Cd(x,aO)-2/{Q-i)(i + o(l)), 

where C = C{Q). Several researchs on the more general equation 

AU = p{x)f{U) 

have been done with different assumptions on / and on the weight p, with asymptotic expansions 
near dn , see for instance [2],[3],[7],[9],[l6],[l7],[l9],[2D],l22];see also [T], [IQ] for quasilinear 
equations. These results rely essentially on the comparison principle valid for this equation, 
and the construction of supersolutions and subsolutions. 

The existence and the behavior of solutions of (|1.6|) in Q\ {0} which blow up at 0: 

lim U (x) = oo, 

called large (or singular) at 0, have also been widely investigated during the last decades, see 
for example [23], and the references therein. There exists a particular solution in ]R^\ {0} 
whenever Q < N/{N - 2) or iV = 1, 2, given by U*{x) = C* jxp^/^^"^) , with C* = C*(Q, N). 
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If Q > N/{N — 2), there is no large solution at 0, and the singularity is removable. If Q < 
N/{N - 2) or iV = 2, any large solution satisfies lim|3,|_^o \xf^^'^~^^ U = C*, or 

lim |a;|^"^ [/ = a > if iV > 2, lim jln |x| | C/ = a > 0, if iV = 2. (1.7) 

|a;|->0 |x|^0 

There exist solutions of each type, distinct from U*. Moreover, up to a scaling, there exists a 
unique positive radial solution in K'^\ {0}, such that ()1.7p holds and limu|__^oQ Ixl ' ^ ' U = C*, 
see ^23j and also [1]. 

In Section [2] we consider the blow up problem of system (jl.ip at the boundary. 

Up to our knowledge all the known results for systems are related with systems for which 
some comparison properties hold, for example 

An = u'^v^, 
Av = uf'v"', 

where s,m > 1, 5, fj, > 0, and Sfj, < {s — l){m — 1), of competitive type, see [13], or 5, < 0, 
of cooperative type, see [8]; see also some extensions to problems with weights in [21], or with 
quasilinear operators in [14], [24], |25j . and cooperative systems of Lotka-Volterra in |12j . 

On the contrary the problem ()l.ip - (jl.2p has been the object of very few works, because it 
brings many difficulties. The main one is the lack of a comparison principle for the system. 
As a consequence all the methods of supersolutions, subsolutions and comparison, valid for the 
case of a single equation fail. 

Until now the existence of large solutions is an open question in the nonradial case. In the 
radial case the problem was studied in [i5j, without weights: a = b = 0. It was shown that 
there are infinitely many nonnegative radial solutions to (jl.ip which blow up at the boundary 
of a ball provided that ()1.4p holds, and no blow up occurs otherwise. In particular, there exist 
solutions even in the case where either u or w vanishes at 0. This shows the lack of a Harnack 
inequality, even in the radial case. The precise behavior of the solutions was obtained in 
forA^ = l,a = 6 = 0, where system (jl.ip is autonomous, with an elaborate proof wich could 
not be extended to higher dimension. 

Our first main result solves this question in any dimension, with possible weights, and 
moreover we give an expansion of order 1 of the solutions: 

Theorem 1.1 Let {u,v) be any radial nonnegative solution of defined for r E (ro,R), 

ro > 0, unbounded at r = R. Then lim u(r) = lim v{r) = oo, and u,v admit the following 

r^R r^R 

expansions near R : 

u{r) = Aid{r)-^{l + 0{d{r))), v{r) = Bid{r)-^1 + 0{d{r))), (1.8) 
where d{r) = R — r is the distance to the boundary, and 

Ai = (7(7 + l)m + 1))')'/'', B^ = m + 1)(7(7 + 1))^)^/^. (1.10) 
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Our proof is essentially based on a new dynamical approach of system (jl.ip . initiated in [3]: 
we reduce the problem to a quadratic, in general nonautonomous, system of order 4, which, 
under the assumptions of Theorem 11.11 can be reduced to a nonautonomous perturbation of 
a quadratic system of order 2. We then show the convergence of the solution of the original 
system to a suitable fixed point by using the perturbation arguments of [TH] . 

Theorem 1 1 . 1 1 can be applied to sign changing solutions of some elliptic systems, in particular 
to the biharmonic equation, where 6=1: 

Corollary 1.2 Let /i > 1, 6 G M. Then any radial solution u of the problem 

A'^u = \x\^ lulf" in (ro,i?), u{R) = oo, (1.11) 

satisfies 

u{r) = A(i(r)-^/(^-i)(l + 0{d{r))), (1.12) 
with A^'-'^ = 8(/i + 3)(^ + l)(3/i - 1)(m - 1)"^. 

We notice here a case where we find an explicit solution: for > 4 and jj, = equation 
A^n = admits the solution in the ball -6(0, 1), 

andv = Au = C{N -4:){l-r^)-^/^{N -2r^) > 0, and ^B) and (frT2]l hold with 7 = = 

TV 
2 • 

In Section [3] we consider the problem of large solutions at the origin, that is (jl.ip - ()1.3p . 

System (II. ip admits a particular radial positive solution {u*,v*), given by 

u*(r) = A^vr"'^"''', v*{r) = Bnv'^''-'', r = \x\, (1.13) 

where 

(2 + a) + (2 + 6)^ ^ (2 + 6) + (2 + a)/. 

7a,6 = ^ > 0, Ca,b = > 0, (1.14) 

= -fa,b{la,b -N + 2) {Ca,b{Ca,b -N + 2))^ , B§ = Ca,b{Ca,b - N + 2) (7a,fe(7a,6 -N + 2)^ , 

whenever 

mm{ja,b,L,b}>N-2, or TV =1,2. (1.15) 
Note that in particular 70,0 = 7, £.0,0 = 

The problem has been initiated in [26j and [S], see also |27j. Let us recall an important 
result of [5] giving upper estimates for system (jl.ip in the nonradial case, stated for N > 3, 
but its proof is valid for any > 1. It is not based on supersolutions, but on estimates of the 
mean value of u, f on spheres: 

Keller- Osserman type estimates [5]. Let Q be a domain of {N > 1), containing 0, and 
u,v € C^(r2\{0}) be any nonnegative subsolutions of lll.l]) . that is, 

-Au+ |x[V < 0, 
-Av + < 0, 
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with iJ,,6 satisfying (L^. Then there exists C = C{a,b, 5, fi, N) such that near x = 0, 

-(/(x) < C ^(x) < Cla;p^"•^ (1.16) 

Moreover, one finds in |5j a quite exhaustive study about all the possible behaviors of the 
solutions (radial or not) in r2\{0}. 

Here we complete those results by proving the existence of local radial solutions large at 
of each of the types described in [5], see Propositions 13.21 13.41 in Section [3l By using these 
results, we obtain our second main result in this work, which is the following global existence 
theorem: 

Theorem 1.3 Assume that N > 2 and that I11.15\) holds. Then there exists a radial positive 
global solution of system in ]R^\{0}, large near 0, unique up to a scaling, such that 

lim r'^'^.by^ ^ Yi-^ ^(,a,by ^ j^^. (1-17) 



r— >oo r— >oo 



and, for N > 2, and up to a change of u, /x, a, into v, 6, b, when 6 < it satisfii 



lim "^u = a > 0, 



es 

limr^-^v = /3>0, «/^<f^, 
limr(^-2tM-(2+;>)^ = /3>0, iffi>^, 
limr^-2|lnrr% = (3 > 0, if = §±^, 



and for N = 2, 

lim llnrp"*^ li = a > 0, lim llnrp"*^ v = (3 > Q. 

r— >-0 r— i>0 

Our proof also relies on the dynamical approach of system (jl.ip in dimension A'^ by a 
quadratic autonomous system of order 4, given in [4J. Finally we give an application to the 
biharmonic equation: 

Corollary 1.4 Let N > 2. Assume that 1 < fi < ^jy,^^^ • There exists a positive global 
solution, unique up to a scaling, of equation 

A\ = |x|V 

in {0} , such that 

lim = a > 0, lim / (''-^^ = C, 

r— >0 r— >oo 

where C'^ = {A + b){N + 2 + b - [N - 2)^) (2^ + 2 + 6)(7V + 6 - (TV - 4)^) (// - 1)-^. 

2 Large solutions at the boundary 

This section is devoted to the study of the boundary blow up problem for nonnegative radial 
solutions of (jl.ip . We begin by observing that system (jl.ip admits a scaling invariance: if (n, v) 
is a solution, then for any 9 > 0, 

r ^ {e^'^-bu{dr),e^^'^v{dr)), (2.1) 
where 7a,6,Ca,fe are defined in (jl.l4p . is also a solution. 
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2.1 Existence and estimates of large solutions 

We say that a nonnegative solution (u, v) of (jl.ip defined in (0, R) is regular at if v G 
{0,R) n C{[0,R)). Then u,v e C^{[0,R)) when a,b > -1, and moreover u'{0) = v'{0) = 
when a,b > —1, and n, w G C^([0, i?)) when a,b > 0. 

We first give an existence and uniqueness result for regular solutions: 

Proposition 2.1 Assume il.5\) and only that D = 5fi — 1 0. Then for any uq, vq > 0, there 
exists a unique local regular solution {u,v) with initial data {uq,vq). 

The result follows from classical fixed point theorem when uq,vq > 0, by writing the problem 
in an integral form: 

u{r) =uo+ r r^-^ r e^-^+%\e)de, v{r) = vo + f r^"^ f 0^-^+V(0)ci0. 



In the case uq > Q = vq, the existence can be obtained from the Schauder fixed point theorem, 
and the uniqueness by using monotonicity arguments as in [15]. We give an alternative proof 
in Section [31 using the dynamical system approach introduced in [1] , which can be extended to 
more general operators. 

Next we show that all the nontrivial regular solutions blow up at some finite i? > 0, and 
give the first upper estimates for any large solution. Our proofs are a direct consequence of 
estimates (|1.16p . 

Proposition 2.2 (i) Assume {1-4^ and ( fi.5|) . For any regular nonnegative solution {u,v) ^ 
(0,0), there exists R such that u and v are unbounded near R. 

(a) Any solution {u,v) which is nonnegative in an interval {rQ,R) and unbounded at R, 
satisfies 

lim u = lim v = lim u = lim v' = oo. (2-2) 

and there exists C = C{N, 6, fi) > such that near r = R, 

u{r) <C{R-ry^, v{r) < C{R - r)'^ . (2.3) 

Proof, (i) Let (u, v) be any nontrivial regular solution. Suppose first that vq > 0. Then from 
p.ip . r^~'^u' is positive for small r, and nondecr easing, hence u is increasing. If the solution is 
entire, then it satisfies (ll.lGh near oo: indeed by the Kelvin transform, the functions 

u{x) = \x\ u[x/ \x\ ), v[x) = \x\ v[x/ \x\ ), 

satisfy in -6(0, 1)\ {0} the system 

j -Au+ \xfv^ = 0, 
I -Av + \xfu^' = 0, 
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where a = {N - 2)6 - {N + 2 + a),b = {N - 2)fi - {N + 2 + b), and jaM, Ca,b are replaced by 
N — 2 — ^a,b, N — 2 — ^a,b- Then the estimate (jl.l6p for (n, v) imphes the one for (u, v) and thus 
u tends to at oo, which is contradictory. Furthermore, from 

"-"°+ (2 + aK"lr + a) "'' " ^ ^° + (2TWT6)^'' 
u and V blow up at the same point > 0. 

(ii) Since r^~^u' is increasing, it has a limit as r — )• i?. If this limit is finite, then u' is 



bounded, implying that u has a finite limit; this contradicts our assumption. Thus (j2.2p holds. 
By ()2.ip we can assume i? = 1 and make the transformation 

\ e-^ if TV = 2, ^ ^ 

(in particular r = 1 — s if = 1), so that s describes an interval (0, sq], sq > 0, and we get the 
system 

Uss = F{s)v^ , . 

vss = G{s)u^ ^^-^^ 

with 

F{s) = r^^-^+-, G{s) = r^^-^+b- (2.6) 



hence lims_>.o F = limg^Q G = 1. Then 



\ -Vss + 



Vss + < 

in some interval (0,si], thus from the Keller-Osserman estimates (jl.l6p . there exists C = 
C{N,5,ix) > such that u{s) < Cs-^,v{s) < Cs'^, near s = and (USD follows. ■ 

2.2 The precise behavior near the boundary 

In this section we prove Theorem II. 1[ 

2.2.1 Scheme of the proof 

Consider a solution blowing up at i? = 1. In the case of dimension = 1, and a = b = 0, we 
have that F = G = 1 in (|2.5p . and we are concerned with the system 

^ ""1 (2.7) 

Vss = W^- 

Following the ideas of [1], we are led to make the substitution 

x(t) = -^, y(0 = -^, z(t) = < wit) = '-^, 

U V Us Vs 
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where t = \ns, t describes (— cx),to]) ^-nd we obtain the autonomous system 



f2i 



( Xt = X[X + l + Z], 

Yt = Y[Y + l + W], 

Zt = Z[l-5Y-Z], 

^Wt = W[l- iiX -W]. 

We study the solutions in the region where X, y > and Z,W <0. In this region system ()2.8 
admits two fixed points 



O = (0, 0, 0, 0), Mo,i = (7, -1 - 7, -1 - 



(2S 



where 7 and ^ are defined in ()1.9p . We intend to show that trajectories associated to the large 
solutions converge to Mo,i. Observe that system (j2.7p has a first integral, which is a crucial 
point in what follows: 



,<5+l 



UsVs 



fi + l 6 + 1 



C, 



equivalently 



2i XZ YW . 

+ 1 fl + l 



C. 



Since any large solution at r = 1 satisfies limr-^i u = limr_>i v = 00, we obtain 



as t 



5 + 1 fi + l 

-00. Thus, eliminating W, we get the nonautonomous system of order 3 

X[X + l + Z], 

Y[Y + l]-{^l + l)XiY + ^,) + o{e^'), 
Z[1-6Y - Z]. 

which appears as a perturbation of system 




(2.10) 



Yt 
Zt 



X[X + 1 + Z], 

y[y + i]-(M + i)x(y + ^), 

Z[1-6Y -Z]. 



(2.11) 



Moreover, by using a suitable change of variables, system (j2.1Up reduces to a nonautonomous 
system of order 2, and we can show that the last system behaves like an autonomous one. Then 
we come back to the initial system and deduce the convergence. 

In the case > 1 or a, b not necessarily equal to 0, we first reduce the problem to a system 
similar to (12. Sp . but nonautonomous, and we prove that it is a perturbation of (12. 8[) . Moreover 
we produce an identity that plays the role of a first integral, allowing us to reduce to a double 
perturbation of (12. lip . We manage with the two perturbations in order to conclude. 
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2.2.2 Steps of the proof 

Our proof relies strongly in a result due to Logemann and Ryan, see [TS]. We state it below for 
the convenience of the reader. 

Theorem 2.3 JTSi Corollary 4-lJ Let h : R+ x R^^ — > M^-'^ be of Caratheodory class. Assume 
that there exists a locally Lipschitz continuous function h* : M^^ — > M.^'^ such that for all compact 
C C and all e > 0, there exists T > such that 

sup ess sup \ \h{T,x) — h*{x)\\ < e 

c&C T>T 

Assume that x is a hounded solution of equation Xr = h{T,x) on M+ such that x{0) = xq. 
Then the co-limit set of x is non empty, compact and connected, and invariant under the flow 
generated by h* . 

The proof of Theorem [TTT] requires some important lemmas. By scaling we still assume that 
R = l. 

Lemma 2.4 Let {u,v) be any fixed solution of system lil.l]) in [rQ,l), unbounded at 1. Let us 
set t = logs, where s = ^'"^(r) is defined in \2.4\) - Let F,G be defined by h2. 6]) . Then the 
functions 

X{t) = > 0, Y[t) = > 0, Z{t) = < 0, Wit) = '-^^ < 0, (2.12) 

U V Us Vs 

satisfy the (in general nonautonomous) system 



{ Xt = X[X + l + Z], 

Yt = Y[Y + 1 + W], 

Zt = Z[1-6Y - Z -a{t)], 

Wt = W[l- ^iX -w - p{t)], 



(2.13) 



where 



2N-2 + a . , 2A^-2 + 6 , 

a(t) = —e\ Mt) = —eK (2.14) 

^' l + (iV-2)e* l + (iV-2)e* ^ ' 

Moreover we recover u,v by the relations 

u = s'''F-i}G-i{\Z\X)T}{\W\Y)-n, v = s^^F'DG-Ty{\W\Y)i{\Z\X)D. (2.15) 

Proof. Since {u,v) is unbounded, (j2.2p holds. We make the substitution (j2.4p . which leads 
to system (|2.5p , with F, G given by (|2.6p . Clearly we can assume that Ug < and < on 
(0, So], lims_>.o \us\ = limsH-o \vs\ = linis^o'u = lims^o^^ = oo. Then we can define X,Y,Z,W 
by (|2.12p and we obtain system (|2.13p with 

ait) = -s^, m = 

then (|2.14p follows, and we deduce (|2.15p by straight computation. ■ 

Next we prove that system (j2.13p is a perturbation of the corresponding autonomous system 
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Lemma 2.5 Let N >1. Under the assumptions of Theorem M.li there exist k > and t < to 
such that 

1/k < X,Y,\Z\,\W\ <k for t <t. (2.16) 



Moreover, setting 



,,,, XZ YW voit) , 

(3 + 1 /X + 1 /X+1 



we have w{t) = 0(e*) as t ^ — oo. 

Proof. We establish some integral inequalities, playing the role of a first integral, then we use 
them to prove (|2.16p . and finally we deduce the behavior of vj. 

(i) Integral inequalities. Let cr, € M and set 



HaAs) = r^-^ ( UsVs - F{s)—— - G{s 



r^-^uve-^' 



6 + 1 ' ' + 1 1 + {N -2)s 
( X{Z + a{t)) Y{W + m) \ 

6+1 + ^+1 j' 



where 



a{6 + l)s , , 61(^ + 1)5 

m = , and m- 



l + {N-2)s ^' l + {N-2)s 

It can be easily verified that 

H'^eis) = {N-2-a-e)usVs+F{s)-—{N+a-a{6+l)) + G{s)—-{N+b-e{fi+l)). (2.18) 

+ 1 /U + 1 

By choosing first the constants a = ai > and 9 = 9i > large enough, we obtain that 
H'^^0^{s) < and thus H„^q^{s) > —Ci for some t Ci > 0; next choosing o" = cr2 < 
and 9 = 62 < Q and large enough in absolute value, we obtain that H'^^ $2^^^ ^ ^"^^ thns 
Hcj^fiii^) ^ ^"2 for some C2 > 0. Hence, there exists functions ai{t), (3i{t),i = 1,2, which are 
0(e*) as t ^ —00, and such that 

6 + 1 ^ + 1 uv 

^y^E^_tM}l^m±M>l<c,r-~^'-. (2.20) 

+ 1 /U + 1 uv 

(ii) Estimates from below in 112. 16\) . Using that Ugs < and multiplying by 2us < we obtain 

(ul), > 2v^Us = i2v^u), - 26v^~^uvs > i2v^u), 

since < in (0, sq], hence — 2v^u <C = (u^ — 2v^u){so); since lims^.o v^u = 00, it follows 
that ul < {5/2)v^ u on (0, si] , for sufficiently small si. Using the same method for the second 
equation, we obtain from (I2.12P that 

X{t) <3\Z{t)\, Y{t) <3\W{t)\, on (-cx),ti]. (2.21) 
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Also, from the generalized L'Hopital's rule, 



— \Z\ —sv^ —Ssv^-^Vs + v^ —/SY-1 
lim = lim < lim = lim 



s^O F s~^0—U,^ s^O — U,,s s^O V F 

and by symmetry 

(5Tmry > 1 + Ih^r|Z|, uAh^X > 1 +^\W\. (2.22) 

s^O s^O s-s>0 s->0 

Suppose now that lim^_^_^X = 0. From (I2.22p . limf_j._oo-^ ^ ^/l-i-, hence there is a sequence 
{in} —oo of local minima of X such that lim„_j.oo ^(^n) = 0, and from the definition of 
X in (I2.12p . X{tn) > for all n sufficiently large. At each t„ we have that Xt{tn) = and 
Xtt{tn) > 0. From (|2.13p . using that X{tn) 7^ 0, we have that X{tn) + 1 = |-2^(in)| and hence 
|Z(t„)| > 1. Since Xu{tn) = X{tn)Zt{tn), it follows that Zt(tn) > 0, and thus, from the third 
equation in ([233]) . 1 - (5y(t„) + |Z(t„)| < 0, implying 

l<l + \Z{tn)\<SYitn) + a{tn). (2.23) 
From (12191) and ([221]), we deduce 

< SVmt) + (/^ + 1)^) + 3(/i + 1)^^ + 0(e'*), 

+ 1 

hence limn^ooY{tn) = 0, which contradicts (|2.23p . We conclude that lim ^^_^X > 0, and 
similarly for Y, thus X, Y, \Z\, \W\ are bounded from below. 

(iii) Estimates from above. From (j2.3p . s'^u and s^v are bounded as s — )• 0, thus from (j2.15p and 
IK2T\i . is bounded as t —7- oo. Since X, Y are bounded from below, they are bounded 
from above, and then also \Z\ and \W\, from (j2.22p . hence ()2.16p holds. 

(iv) Conclusion. From (I2l9l) . ([220]), since X,Y are bounded and \ai{t)\, \/3i{t)\ < Ce\ 

^y>^ + ^-C3e' and XY<^ + ^ + C,e\ (2.24) 
+ 1 /i + 1 o + l/u + l 

for some 6*3,(74 > 0. Then we deduce ()2.17p . ■ 

Next we show that a convenient combination of our solution {X, Y, Z, W) satisfies a system 
of order 2. We have 

Lemma 2.6 Under the assumptions of Theorem \l.l[ and with the above notations, let 

= y^-m, r = -[z(.ya. (2.25) 

Then (x, y) lies in the region 

for T > f > 0, and satisfies 



Xr = x{—x — 6y + 2)+wi{T) 
Vr = (5^-y)((5 + l)y-(/i + l)x) + ro2(T), 

where wi{t) = 0{e~^'^) and '0J2{t) = 0(6"^"^) for some K > 0, as t ^ oo. 



(2.26) 



11 



Proof. We first reduce system (|2.13p to a system of order 3: from relation ()2.17p we eliminate 
W in the system (j2.13p and obtain 

Xt = X[X + 1 + Z], 

Yt = Y[Y + l]-{f, + l)X(Y + ^)+wit), 
Zt = Z[1-6Y - Z -a{t)], 

which is a perturbation of system (j2.1ip . Next, defining x = y = we get the system 



Xt = Z 
yt = Z 



x{2 — X — 5y) + zui] 
^-y){{6 + l)y-{fi + l)x)+w, 



with 



VJi = = )' ^2 = ^2 = )' <^2.27) 

from Lemma 12.51 and then 

Zt = Z{1 + Z{5y-l) + a{t)). (2.28) 

and r(t) defined by (j2.25p for t < t describes [0, oo) as t describes (— oo, t], and t /2k < \t\ < 2kT 

for t < i. Hence we deduce (j2.26p . and the estimates of wi,zu2. Notice that < x,y < k'^ 
for any r > from (|2.16p . and from (j2.17p . for r > f > 0, 

1 1 , YW ,,,, ^ 1 

y-^ = w(t^ + «(i))> 



6 + 1 + 1 ^ - 2{fj, + l)k^' 

ending the proof. ■ 

Hence system (j2.26p appears as an exponential perturbation of an autonomous system that 
we study now: 

Lemma 2.7 Consider the system 

= ^(2-x-6y) ^2.29) 



The fixed points of system i2. 29\) are O = (0,0), and 

■ M p (5 + 2 I \ ( 2{5 + l) 2(/i + l) 



S + iJ^ \S + V6 + lJ' ^^'^'^^ \^s + 26 + l' fi6 + 26 + l 

and mo is a sink. Any solution of the system ^2. 29\) which stays in the region TZq converges to 
the fixed point mo as r ^ oo. 

Proof. The point mo is a sink: the eigenvalues of the linearized system of ()2.29p at mo are the 
roots £i , £2 of equation 

2 <5/i + 3 + 2/i + 2^ nS + 2^i+J^ ^ 
^ ^i6 + 25 + l ^%i5 + 25 + l ' 

equivalently 

(7 + l)^ + {-i + i + l)l + 2(e + 1) = 0, (2.30) 
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and they have negative real part. Next we show that (j2.29p has no hmit cycle in (0,oo) x 
(l/(<5 + 1), oo). Let B = xP{y - ^)~^ where p, q are parameters. Writing (I2.29[) under the 
form Xr = T{x,y), yt = G{x,y), we obtain 

V • (BiT, g)) = B^xr + B^nVr + BiT^ + Gy) := MB, where 
M =ip-l-p-q{fi+ l))x - ip5 - qi5 + 1) + 35 + 2) (y - -i-) + + + Qi^ + 1) + 1 



(5 + r 5 + 1 



Choosing q = ^^l23+i ^^"^ P = ^^ ~ ^ ~ Qif^ +1)) find that 



Hence, by the Bendixson-Dulac Theorem, system (j2.29p has no limit cycle. From the Poincare- 
Bendixon Theorem, the w-limit set F of any solution of (j2.29p lying in TZq is fixed point, of a 
union of fixed points and connecting orbits. But is the unique fixed point in TZq. Then any 
solution in TZq converges to uiq as r ^ oo. ■ 

Remark 2.8 It is easy to prove that there exists a connecting orbit joining the two points io 
and rn-o, but it is not located in TZq. 

We can now conclude. 

Proof of Theorem I I.IL (i) Convergence for system /i2.26\) . From Proposition 12.31 the w-limit 
set T, of our solution {x,y) of ()2.26p is nonempty, compact, connected and contained in TZq, and 
S = U^gE T>f ■^)' where f{T,i) denotes the trajectory of (I2.29j) such that ^p{f,i) = I. Since 
limT-_!.oo f{T,i) = ruQ, there holds rriQ € S. Since niQ is a sink of (j2.29p . then from the standard 
stability theory, see for example ^ Theorem 3.1, page 327], {x,y) converges to mo. 

(ii) Convergence for system i2.1&]) . By setting g(t) = 1/Z, we find from (|2.28p that g'+{l—a)g = 
1 — 6y, hence by L'Hopital's rule, 

z= lim '^f -">' = lim = -^^^±^ = + , ) = Zo. 

t->--oo t^-oo fg^Jti^^'^^y t^-oo I — Oy 110 — 1 

Hence 

2(5 + 1) 2(/i + l) ^ 

lim X = — lim xZ = — = j = Xq, lim Y = lim yZ = — = ? = 

t— oo t^-— oo ^0 — 1 t— oo t— oo ^0 — 1 

Finally, from (I2l7l> . we obtain limt^_oo = -(1 + = Wq. That means {X,Y,Z,W) 
converges to Mq^i defined at (j2.9p . Then from (j2.15p we deduce the estimates 

u{r) = Aid-^il + o(l)), v{r) = Bid-^l + o(l)) 

where Ai , Bi are given by and (jl.lOp . 

(iii) Expansion of u and v. We first consider system (j2.26p . Setting x = xq + x,yo + y, we find 
a system of the form 

(xr^yr) = + Q{x,y) + (^1,1172) 
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where (x, y) — > (0, 0) , the eigenvalues £i, £2 of A satisfy max(Re(£i, £2)) = — < —1/(7+1), and 
Q is quadratic and ■cui{t),W2{t) = 0{e~^'^). There exists an eucUdian structure with a scalar 
product where (^(x, y), (x, y)) < — m||(x,y)||^ . Then the function r i-> ^^(r) = ||(x,y)||(r) 
satisfies an inequality of the type rj^ < — (m — e)ri + Ce~^'^ for any e > and r large enough. 
Then 

r,{T) = 0{e-^^) + 0(e-(™-")"). (2.31) 

Then the convergence of (x,y) to (xo,2/o) is exponential. From (|2.28|) . the convergence of Z to 
Zq is exponential. Writing r under the form 

{Zo - Z), 

we deduce that r = c + ZqI + 0{e^^) for some A: > 0. From (j2.27p we obtain that 'dji,W2 = 
O(e-^o^) with Kq = 1/\Zq\ ; taking K = Kq = 1/(7 + 1) in (fOTT) . we find that ?7(r) = 
O(e-^o^) = 0(e*), because m > ETq- Then from (l2:25]l we deduce that \Z - Zq\ = 0(e*), and 
then from dOS]), \X - Xq\ + \Y - Ya\ = 0(e*), and in turn \W - Wq\ = 0(e*) from (I2T7D . 
Finally we come back to u and v by means of (j2.15p : recalling that s = e* and r = 1 + 0(s) as 
s ^ 0, we deduce that 

u{r) = ^is-^(l + 0(s)), v{r) = Bis'^il + 0{s)) 

and the expansion (jl.Sp follows from (j2.4p . ■ 

Proof of Corollary II. 2L Let u be a radial solution of (jl.lip . Then u and v = satisfy 

{Au = V 
Av = 

and then u{r) > in (ro,i?) and u{R) = 00. Integrating twice the second equation in this 
system, we have that lim v{r) = 00 and Theorem 1 1 . 1 1 applies . ■ 



2.3 The set of initial data for blow up 

Here we suppose a = 6 = 0. By scaling, for any p > there exists solutions which blow up at 
p. Let us call p{uq,vq) the blow-up radius of a regular solution with initial data {uq,vq). From 
(|^ . we find 

p(A'^uo, A^wo) = ^'^p{uq,vq). 

Then for any {uq,vq) € there is a unique A such that p{\"'uo,\^vq) = 1. Thus there exist 
infinitely many solutions blowing up at i? = 1, including in particular two unique solutions with 
respective initial data (no,0) and (0, uq). Using monotonicity properties, it was shown in [15] 
that the set 

5 = s (uo, vq) € [0, 00) X [0, 00) : lim u = lim v = oo\ 

is contained in [0, uq\ x [0, vq]. Next we give some properties of S extending some results of [15] 
to higher dimensions. 

Proposition 2.9 Let N >l. Ifuim{5,p} > 1, then S is a simple curve joining the two points 
(no, 0) and (0, vq). 
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Proof. We claim that the mapping (uq, vq) € [0, oo) x [0, oo) \ {0, 0} i — > p{uo,vo) is continuous. 
As in [111 this will follow from our global estimates. 

(i) The function p is lower semi- continuous. Indeed the local existence is obtained by the 
fixed point theorem of a strict contraction, since min{5, //} > 1, then we have local continuous 
dependence of the initial conditions, even if uq = or vq = 0, and the result follows classically. 

(ii) The function p is upper semi- continuous. We can start from a point tq > instead of 0. 
We prove that for any positive {uo,vo), considering any solution (n, {;) equal to (tto,'L'o) at ro, 
with blow-up point p, for any f > p, any solution (n, v) starting from tq with data sufficiently 
close to {HqjVo), blows up before r : suppose that it is false, then there exists a sequence of 
positive solutions {un,Vn), with data {un,Vn) at tq, tending to {uq^vq), increasing, and blowing 
up at pn > r. We can assume f = 1. Making the change of variables (12. 4p we get solutions of 
system ()2.5p in (0, so], satisfying Co = Co(ro, A^, a, b) 

-Uss + Cqv^ < 
-vss + Con^ < 

with u and v decreasing. In fact estimates (|1.16p hold with a universal constant, in any 
B{0,k)\{0} C such that the mean values of u and v on dB(0,r) are strictly monotone. 
Then there exists a constant C = C{Co,N,5, p) such that 

Unis) < Cs~"' , Vn{s) < C s~^ for s < So, 

that means 

Un{r) < C{r^-^ - Vn{r) < C{r^-^ - l)"? for r e [ro, 1) . 

Passing to the limit we find that u,v are bounded at the point p < 1, which is contradictory. 
Then the claim is proved. Thus 5 is a curve with 

{uq,vq) = p'^ {cos 9 , sin 9) cos 6 , p^ {cos 9 , sin 6) sin 6 , 6's[0,7r/2], 

as a parametric representation. ■ 

3 Behavior of system (11.11) near the origin 
3.1 Formulation as a dynamical system 

In [3] the authors study general quasilinear elliptic systems, and in particular the system 

-Au = -{urr + ^f^Ur) = eirV, , . 

-Av = -{Vrr + ^Vr) = EarV, ^""^^ 

where ei = ±1, £2 = il- Near any point r where u{r) 7^ 0,u'{r) ^ and v{r) ^ 0, v'{r) ^ 0, 
they define 

X(t) = -^, r(t) = -^, Z{t) = -e/-^, w{t) = -e/-^, (3.2) 

U V Ur Vr 
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with t = lnr, so system p.ip becomes 



Xt = X[X-{N-2) + Z], 

Yt = Y[Y -{N -2) + W], 

Zt = Z[N + a-5Y -Z], 

I Wt = W[N + b-fiX-W]. 

One recovers u and v by the formulas 

u=r-^a,t^ \ZX\)^/^{ \WY\Y/^, v=r-^-''' \WY\y/^ \ZXf'^ 

and we notice the relations 7^ f, + 2 + a = (5^^ {, and ^^.6 + 2 + 6 = /U7a f,. 



(3.3) 



(3.4) 



As mentioned in |4J, system (j3.3p is independent of Ej, i = 1,2, and thus it allows to 
study system (|3.1|) in a unified way. In our case ei = £2 = — Ij then XZ = —r°-~^'^v^ /u and 
YW = r^^'^u^ /v, thus we are led to study ()3.3p in the region 

71 = {{X, Y, Z, W) \ XZ < 0, YW < 0}. 

This system is quadratic, and it admits four invariant hyperplanes: X = 0,Y = 0, Z = 
0,W = 0. The trajectories located on these hyperplanes do not correspond to a solution of 
system (j3.ip . and they are called nonadmissible. System p.3p has sixteen fixed points, including 
O = (0, 0, 0, 0). The main one is 



Mo = {Xo,Yo, Zo, Wo) = (7a,fe, ^b, N-2- 7^,^, N 



which is interior to TZ whenever (jl.lSp holds; it corresponds to the particular solution {u*,v*) 
given in (I1.13p . Among the other fixed points, as we see below. 

No = {0,0,N + a,N + b), 

Ro = {0,-{2 + b),N + a + {2 + b)S,N + b), So = (-(2 + a), 0, iV + a, iV + 6 + (2 + a)^) , 
are linked to the regular solutions, and 

^0 = (A^-2,iV-2,0,0), Go = (N -2,0,0, N+b-{N-2)fi), Hq = {0, N -2, N+a-{N -2)5,0), 



Po = {N- 2, {N-2)fi-2-b,0,{N + b-{N- 2)fi)), 
Qq = ([N -2)6-2-a,N -2,N + a-{N - 2)6, 0), 

and Moare linked to the large solutions near 0. Notice that Pq ^ TZ for ^ < fi < §^ and 
Qo 7^ for < 6 < We are not concerned by the other fixed points 

lo = (N -2,0,0,0), Jo= iO,N -2,0,0), Ko= (0,0, N + a,0), Lo= (0, 0, 0, iV + 6), 

which correspond to non admissible solutions, from and 

Co = {0,-{2 + b),0,N + b). Do = {-{2 + a),0,N + a,0), 

which can be shown as non admissible as t — >■ —00. 
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3.2 Regular solutions 

First we give an alternative proof of Proposition 12. 1[ 

Proposition 3.1 Assume il.5\) and D ^ 0. Then a solution {u,v) is regular with initial data 
{uq,vq),uo,vo > (resp. {uo,0),uo > 0, resp. {0,vo),vo > 0), if and only the corresponding 
solution {X,Y, Z,W) converges to Nq (resp. Rq, resp. So) as t ^ —oo. For any uq,vo > 0, not 
both 0, there exists a unique local regular solution {u,v) with initial data (uqjVq). 

Proof. The proof in the case uo,vq > is done in jl^ Proposition 4.4]. Suppose uq > = vq, 
and consider any regular solution {u,v) with initial data {uo,0). We find 



r'+\l + o{l)), 



N + b 



{N + b){2 + b) 



r2+\l + o(l)), 



„.fJ-\N-l+a+{2+b)5 



[r'^^-'u')' 



{{N + b){2 + b)Y 



(1 + 0(1)), u' 



„ 5M„l+a+(2+fe)5 

Uq / 



{{N + 6)(2 + b)Y{N + a + (2 + b)5) 



then from (|3.2p the corresponding trajectory {X, Y, Z, W) converges to Rq as t 
we show that there exists a unique trajectory converging to Rq. We write 



-oo. Next 



Ro 



{0,Y,Z,W) = {0,-{2 + b),N + a + {2 + b)6,N + b). 



Under our assumptions it lies in K. Setting Y = Y + Y,Z = Z + Z,W = W + W,t\ie 
linearization at Rq gives 



Y 



Y + W 



Xt = AiX, Yt 
the eigenvalues are 

Ai = 2 + a + (5(2 + 6) >0, A2 



Zt = Z 



-5Y - Z 



Wf = W 



(2 + 6) < 0, X3 = -Z < 0, A4 



-fiX -W 



-{N + b) < 0. 



The unstable manifold Vu has dimension 1 and Vun{X = O} = 0, hence there exist precisely 
one admissible trajectory such that X < and Z > 0. Moreover it satisfies 

lim e"^i*X = Ci > 0, lim Y = Y, lim Z = Z, lim W = W. 

t— 00 t— >— 00 t— )■— 00 t— )■— 00 

Then from (13. 4j) u has a positive limit uq, and v=0{e^^), thus v tends to 0; then {u, v) is regular 
with initial data (uq , 0) . By (j2.ip we obtain existence for any [uq , 0) and the uniqueness still 
holds. Similarly the solutions with initial data {0,vq) correspond to Sq. ■ 



3.3 Local existence of large solutions near 

Next we prove the existence of different types of local solutions large at 0, by linearization 
around the fixed points Aq,Gq, Hq, Pq,Qq. For simplicity we do not consider the limit cases, 
where one of the eigenvalues of the linearization is 0, corresponding to behaviors of u, v of 
logarithmic type. All the following results extend by symmetry, after exchanging u, 6, a, ^a,b 
and v,fj,,b,S,a,b- 



17 



Proposition 3.2 Assume N > 2. 

(i) If 5 < and fi < then there exist solutions {u,v) to lil.l]) such that 

= /3 > 0. 



lim = a > 0, 



lim r 

r-i>0 



(3.5) 



If ^ > or fJ. > there exist no such solutions. 



(a) Let ja,b > N — 2 and let fi < or > Then there exist solutions {u, v) of ( li. jj) 

such that 



lim r 



N-2 



u = a>0, lim r(^-^)f^-i^+b)y = Ma) > 0, 



(3.6) 



with I3{a) = ai'/iiN - 2)n - N - b){{N - 2)n - 2 - b). If ja,b < N -2, there exist no such 
solutions. 



(Hi) If fi < then there exist solutions {u,v) of such that 



Ifli> 



2+b 
N~2 



lim r^ ^li = a > 0, 



there exist no such solutions. 



lim V = B > 0. 



(3.7) 



Proof, (i) We study the behaviour of the solutions of (j3.3p near Aq as t ^ —oo. The lineariza- 
tion at ^0 gives, with X = N-2 + X,Y = N-2 + Y, 



X + Z 



Yt = iN- 2) 



Y + W 



Zt = X3Z, Wt = X4W, 



Xt = {N- 2) 
with eigenvalues 

Ai = A2 = (iV - 2) > 0, A3 = iV + a - (iV - 2)5, = N + b - {N - 2)fi. 



Ii6 < 



N+a 
N-2 



and fi < 



N+b 
N-2' 



then we have A3, A4 > 0; the unstable manifold Vu has dimension 4, 

-cxj, interior to TZ, then 



then there exists an infinity of trajectories converging to Aq as t - 
admissible, with Z,W < 0. The solutions satisfy lim e~^''^Z = Zq < and lim e~'^'**VF = 

t—^—oo f— 00 

Wq < 0, with lim X = lim Y = N — 2. Hence from ()3.4p . the corresponding solutions (n, v) 

i— 00 t— 00 

of ([H]) satisfy If 6 > or > then A3 < or A4 < 0, respectively, and V„ 

has at most dimension 3, and it satisfies Z = or W = respectively. Therefore there is no 
admissible trajectory converging at —00. 

(ii) Here we study the behaviour near Pq. Setting Pq = {N — 2, y*, 0, VT*), with 
K = (iV - 2)^ - 2 - 6, W, = N + b-{N- 2)fi, 
the hnearization at Pq gives, with X = N-2 + X,Y = Y^ + Y,W = W^ + W, 



Xt = {N- 2) 



X + Z 



Yt = Y, 



Zt = X^Z, Wt = 



-fiX -W 



Y + W 

By direct computation we obtain that the eigenvalues are 

Ai = iV-2>0, A2 = n, X3 = N + a-5Y, = D{^a,b-iN-2)), X^ = -W. 
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Assume first that -jafi > N — 2. Then A3 > 0. If ^ > ^35, then also A2, A4 > and thus Vu has 
dimension 4, then there exist an infinity of admissible trajectories, with Z < 0, converging as 
t — )■ — cxD. If /i < then A2, A4 < 0, thus Vu has dimension 2, and V^nlZ = 0} has dimension 
1, thus there also exist an infinity of admissible trajectories with Z < converging when 

t -00. Then lim e'^^^Z = C3 < 0, lim X = iV - 2, lim Y = and lim W = W*, 

t— >— 00 t— >— 00 t— >■— 00 t— 00 

thus ([331), (-u, satisfy ([M])- If 7a,fe < iV - 2, then A3 < and V„ = V„ n {Z = 0} and there 
is no admissible trajectory converging when t —00. 



(iii) We consider the behaviour near Gq. The linearization at Go gives, with X = N 
X^W = N + b- {N -2)fi + W, 



2 + 



Xt = iN- 2) 



x + z 



Yt = {2 + b-iN- 2)fi)Y, 



Zt = {N + a)Z, Wt = {N + b-{N - 2)fi) 



-fxX -W 



and the eigenvalues are 

Ai = iV-2>0, A2 = 2 + 6-(iV- 2)fi, A3 = iV + a > 0, X4 = (N - 2)n - N - b. 



If /i < 7v^' then A2, A4 < 0. Then Vu has dimension 3, and Vu n {Y = 0} and Vu r) {Z = 0} 
have dimension 2. This implies that Vu must contain admissible trajectories such that X > 



(because iV - 2 > 0), F < 0, Z < and > (because iV + 6 - (A^ - 2);U > 0). Clearly, 
lim A = iV - 2 and Mm W = N + b - {N - 2)^jl > 0. Moreover, hm e'^^iy = C2 < 

t— ^— 00 t— >— 00 t—>-—co 

and hm e-^3*Z = C3 < 0, thus (EZl) follows from ([331) • Let now // > jp^, so that A2 < 0. If 



t—^—ca 
N+b 
Af-2' 



then A4 < 0, Vu has dimension 2, and also Vu n {F = 0}, hence Vu = fl {Y = 0}, 
and there exists no admissible trajectory. If ^ > ^^1) then A4 > 0, Vu has dimension 3 and 
also Vu n {Y = 0}, there is no admissible trajectory. ■ 



Remark 3.3 If fJ, > in (H) the two functions u, v are large near 0. If < j^^, then u 

large near and v tends to 0. 

Next we study the behavior near Mq, which is the most interesting one. 



IS 



Proposition 3.4 Assume N > 1 and I11.15\) . Then (up to a scaling) there exist infinitely many 
solutions defined near r = such that 

r— >-0 r— >-o 

Proof. Setting X = Xq + X ,Y = Yq + Y , Z = Zq + Z ,W = Wq + W , linearized system is 



r Xt = Xo{X + Z), 

Yt = Yo{Y + W), 

Zt = Zq{-5Y-Z), 

Wt = Woi-fiX-W). 
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As described in [Ij, the eigenvalues are the roots Ai, A2, A3, A4, of the characteristic polynomial 

/ - A Xo \ 

= 5\Zo\ iZol-A 

V ^l^oj \Wo\-X / 

= (A - Xo)(A + Zo)(A - yo)(A + Wo) - 6fiXoYoZoWo, (3.8) 

where we recall that Xq, Yq > and Zq, Wq < 0. We write / in the form 

/(A) = A^ + EoX^ + FoA2 + GoA - Ho, 

with 

= Zo-Xo + Wo-Yo, 
= (Zo - Xo)iWo - Yo) - XoZo - YoWo, 
= -YoWoiZo - Xo) - XoZoiWo - Yo), 
= DXoYoZoWo. 

We note that £^0 < 0, Fq > and 2Go = -Eo [YqZo + XoWq] < 0. From we have Hq > 0, 
hence A1A2A3A4 < 0. Hence there exist two real roots A3 < < A4, with 

A4 >max({Xo,yo,|^o|,|VFoj} 

from (j3.8p . and two roots Ai,A2, which may be real or complex. From the form of /(A) in 
(jS.Sp . we also see easily that if the roots Ai, A2 are real, they are positive. Next we claim that 
ReAi > 0. Suppose ReAi = 0. Then /(ilmAi) = 0, then Gg = EqFoGq + E^Hq, and thus, 
dividing by Eo, 

= Gl- EoFoGo + ElHo = ([YoZo + XoWof + 2 [YoZo + XoWo] Fo - AHo) , 

hence [YqZo + XoWo + Fof = Eg + iHo > E^; but 

YoZo + XoWo + Fo = {Xo - Wo){Yo - Zo) G (0, Fo) 

which is a contradiction. Since ReAi is a continuous function of {S,fi), it is sufficient to find a 
value 5) satisfying (jl.l5p for which it is positive. Taking S = fi, the equation in A reduces 
to two equations of order 2: 

/(A) = (A - Xo)\X - \Zo\)^ - d^Xlzl 

= [A^ - (Xo + \Zo\)\ -{5- l)Xo \Zo\\ [A^ - {Xo + |Zo|)A + (1 + 5)Xo \Zo\\ , 

and Xo + \Zo\ > 0, thus the claim is proved. Then Vu has dimension 3 and Vs has dimension 
1. Hence the result follows. ■ 



Remark 3.5 In the case N = 1, two roots are explicit: A3 = — 1, A4 = 2 + 7 + ^, and Ai, A2 
are the roots of equation 

A2-(l + 7 + ^)A + 2(l + 7)(l + e) = 0. (3.9) 
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The 4 roots are real i/ (1 + 7 + — 8(1 + 7)(1 + ^ 0; ^^^^ means 

{6n + 3 + 2fi + 25f - 8{n6 + 26 + l){fi6 + 2/i + 1) > 0, 

which is not true for 5 = fi, but is true for example when 5/fj, is large enough. The roots of 
equation i'J. y\) and the roots of equation 112. !^(J\) relative to the linearization of system 112. 2y\) at 
ruQ are linked by the relations ii = Xi/ \ Zq \ , £2 = X2/ \ Zq \ . Indeed Mq = Mo,i defined at /i2. 9\} 
satisfies relation \2.n\j with = 0, thus (Xq, Yq,Zq) is a fixed point of system i2.11\) and the 



linearization of 12.11\) at this point gives the eigenvalues — 1, Ai, A2. The point is the image 
of {Xq,Yq,Zq) by the transformation i2.25\) . which divides the eigenvalues by \Zq\, due to the 
change in time t ^ t. 



3.4 Global results 

Here we prove our second main result. 

Proof of Theorem 11.31 From the proof of Proposition 13.41 the hnearization at Mq admits 
a unique real eigenvalue A3 < 0. From (j3.8p a generating eigenvector {ui,U2,u^,U4) satisfies 
< and 1*2^4 < 0, and hence it is of the form u = (— a^, — /3^, o"^, p^), or —u. There 
exist precisely two trajectories Tjj and 71^ converging to Mq as t — )• 00 and the convergence 
of X, Y, Z, W is monotone near t = 00; from ()3.4p . the corresponding solutions (u, v) of system 
p^D satisfy (fTTTl) . 

We consider the trajectory 7^ corresponding to u. Let us show that the convergence is 
monotone in all R. Notice that neither of the components can vanish, since system (jl.ip is of 
Kolmogorov type. Near t = 00, X and Y are increasing, and Z,W are decreasing. Suppose 
that there exists a greatest value ti such that X has a minimum local at ti, hence 

Xu{ti) = X{ti)Zt{ti) > 0, = N-2- X{ti), 

thus Zt{ti) > . Then there exists t2 > ti such that Zt{t2) = 0, and 

Zu{t2) = -6Z{t2)Ytit2) < 0, Z{t2) = N + a- 5Y{t2), 

then Yt{t2) < 0. There exists t^ > t2 such that ^^(^3) = 0, and 

Yttits) = Y{ts)Wt{t3) > 0, Y{ts) = N-2- Wits). 

There exists t4 > t^ such that Wt{t4) = and Wuiti) = —W{t4)Xt{t4) < 0. From the definition 
of ti, this implies t^ = ti, and then all the conditions above imply that (X, Y, Z, W){ti) = Mq, 
which is impossible. Hence X stays strictly monotone, and similarly Y, Z, W also stay strictly 
monotone. Since X,Y > 0, and Y,Z < 0, then Tn is bounded, hence defined on M and 
converges to some fixed point L = {li,l2,h,h) of the system as t ^ —00 and necessarily 
h < Xo,l2 < Yo,k > Zo,U > Wo- 

• Case N > 2. First we note that along we always have X,Y > N — 2. Indeed, if at 
some point t we have X(t) = N — 2, then Xt{t) = (N — 2)Z{t) < 0, which is contradictory. 
Hence the possible values for L are Ao, or Pq when ^ > or Qo when 6 > since Iq is 

nonadmissible. By hypothesis, 7a,fe > N — 2, then either n < or S < We can assume 
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that 6 < Then Qo 7^, then L = Aq or Pq. When /i < f±|, then L = Aq. When 

fj, > from Proposition I3.2r i). we have L / thus L = Pq. In the hmit case /i = 

we find Pq = ^o- From the hnearization at Aq we have 

Ai = As = iV - 2 > 0, A3 = TV + a - (iV - 2)(5 > 0, A4 = 0. 

Coming back to the proof of Proposition I3.2l fi). we find that the convergence of Z and X = 
X — (N — 2) to are exponential. From the fourth equation in we see that Wt + > 0, 
hence —1/W < C\t\ near —00. Then, there exists m > such that 

Wt = W^{-1 - fiW-^X) = W^{-1 + 0(e"*); 

integrating over (t,to)) *o < 0, we obtain that W{t) = t^^ + 0(t~^). In turn we estimate Y; 
setting Y = Y + W, then Yt = {N - 2)Y + Y{Y -W) + W{-^iX - W), and thus 

Yt = {{N -2) + e{t))Y + 0{t-^), 

implying Y = Oit""^) and thus Y = N -2-1^'^ + Oit^'^). Next we find that Zt/Z = A3 + + 
0(t~^), which yields l\mt^-oo e~^''^^\t\^^\Z\ = C > 0. Finally, by replacing in (13. 4j) . and deduce 
the behavior of u and v as claimed: 

lim r^~^u = Ci > and lim r^'^] log(r)|-^t; = C2 > 0. 

• Case N = 2. Then necessarily L = O = (0,0,0,0). The eigenvalues of the linearized 
problem at this point are 0, 0, 2 + a, 2 + 6. Since Zt = Z(2 + a — SY — Z) and Y and Z tend to 
as t tends to —00, Z converges exponentially to 0, and similarly W. Since Xt < X'^, it follows 
that X > C \t\^^ near —00. Then 

Xt = X^{1 + Z/X) = + 0(6"^*)) 

for some m > 0, hence X = —1/t + 0(t^^), then the function t 1— )• = u{t)/t satisfies 
(ft/^ = 0(t^^), then if has a finite limit, hence u(r)/lnr has a finite positive limit, and 
similarly for v. ■ 
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